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Abstract 

In this paper we study the wedge disclination within the elastoplastic defect theory. 
Using the stress function method we found exact analytical solutions for all characteristic 
fields of a straight wedge disclination in a cylinder. The elastic stress, elastic strain, elastic 
bend-twist, displacement and rotation have no singularities at the disclination line. We 
found a modified stress function for the wedge disclination. 
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Disclinations are very important lattice defects analogous to dislocations. A disclination is 
characterized by a closure failure of the rotation for a closed circuit round the disclination 
line. There are wedge and twist disclinations. It seems that very strong plastic distortions are 
necessary in order to create disclinations in crystals. If the Frank angle (rotation failure) of the 
disclination is a symmetry angle of lattice, then the disclination is called a perfect disclination. 
Such disclinations have been introduced by Anthony [1] and deWit [2,3]. A wedge disclination 
(rotation axis parallel to the disclination line) may be generated by adding or removing a wedge- 
shaped piece of material and re-welding the material. The smallest value of the Frank vector 
is 7r/2 in a cubic lattice and 7r/3 in a hexagonal lattice. If the Frank angle is not a symmetry 
angle of lattice, the disclination is called partial disclination. They play an important role, e.g., 
in building of twin boundaries (see, e.g., [4]). Disclinations correspond, in general, to Volterra's 
distortions of the second kind (see also [5]). Thus, these defects are of rotational type. They 
are different from the so-called Frank's (spin) disclinations which are (elementary) defects in 
liquid crystals. 

In recent articles [6-9] we have studied screw and edge dislocations in the framework of 
field theory of elastoplasticity. This field theory of elastoplasticity can be considered as a gauge 
theory of defects where the defects cause plasticity. The corresponding gauge fields may be 
identified with the plastic distortion. By the help of this theory the elastic and plastic part of 
the total distortion can be calculated. The total distortion is defined in terms of a displacement 
and consists of the elastic and plastic part. In the case of dislocations (see, e.g., [9]) the elastic 
distortion is continuous even in the dislocation core and the plastic part becomes discontinuous. 
The field equation of the elastic stress in an isotropic medium is an inhomogeneous Helmholtz 
equation (see, e.g., [8]) 

(—2 a \ ° 2 ^fl , , 

1 - K AjCTij = (Tij, K = . (1) 
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Here aij is the classical stress tensor and (i is the shear modulus. The coefficient a\ has the 
dimension of a force and k has the dimension of a reciprocal length. It is important to note that 
Eq. (1) agrees with the field equation for the elastic stress in Eringen's nonlocal elasticity [10-12] 
and in gradient elasticity [13-15] where the factor n is called non-locality parameter or gradient 
coefficient. Using the inverse of Hooke's law for cry and <jjj, Eq. (1) implies an inhomogeneous 
Helmholtz equation for the elastic strain 

(l- K - 2 A)E ij =E ij , (2) 

o 

where E^ is the classical strain tensor. It worth noting that Eq. (2) is analogous to an equation 
for the elastic strain in gradient theory used by Gutkin and Aifantis [13-15]. It is, therefore, 
quite reasonable to use the field theory of elastoplasticity for disclinations. 

Here we consider a straight wedge disclination inside an infinitely long cylinder with outer 
radius R. Disclinations are defined by the disclination line and the Frank vector. The z-axis is 
along the disclination line and coincides with the axis of the cylinder. For a wedge disclination 
the Frank vector is parallel to the disclination line: f2 x = Qj, = and Q. z = 0. 

In absence of body forces, the force equilibrium condition can be identically satisfied by 
using the so-called stress function ansatz [16]. If we specialize to the plane problem, the stress 
function / is related to the stress tensor 

dU -Olyf \ 

-4/ dlJ . (3) 
vAf J 

Here A denotes the two-dimensional Laplacian d 2 x + d yy and v is Poisson's ratio. In addition, 
the strain is given in terms of the stress function as 



/ d 2 yy f - uAf -di y f 



^ [ -% v f dlJ-uAj | . (4) 

We use the classical stress field of a straight wedge disclination in terms of the Airy stress 
function 

dLx -d 2 xy X 

-dlyX dlxX j . (5) 
vAx 

The stress function of a "classical" wedge disclinations is given by 

where r 2 = x 2 + y 2 . That stress function is chosen so that in the case of C = it reproduces 
the stress of a wedge disclination given by deWit [3]. It fulfils the following inhomogeneous 
bipotential (or biharmonic) equation 

AAx = Stt A 5(r). (7) 

More precisely, the stress function x = -Lr 2 lnr is Green's function of the two-dimensional 
bipotential equation. Substituting (3), (5) and (6) into (1) we get 

(A-*)/— ^(tar-j^-c). (8) 
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Now we use the ansatz 

1 - Av 



and obtain 



Its solution reads 



f = ^[^r- w -- ] -C) + f {l) (9) 



(A-4/ a) = -4A(ln,+ ^^-c). (10) 



/a) = ^(i- + ^4^)+W-r«). 



Finally, wc find the solution of (8) 

^ f 2 ( , 1 - 4i/ ^\ , 4 



where the first piece is the "classical" stress function and if„ is the modified Bcsscl function of 
the second kind and n = 0, 1, . . . denotes the order of this function. 

Let us obtain the stress field of the wedge disclination in the cylinder with radius R. It is 
convenient to use cylindrical coordinates 

OVr = - d r f, <7 W = d 2 rr f, U zz = v[p rr + (T w ). (13) 

In this way we find for the non-vanishing components 

= 2^7) { lnr+ 1^2^ + K *W + ^ (2 - ^K^r)) - C), 
o w = { In, + 1 + ^ + K (nr) - ^ ( 2 - «Vtf a (,er)) - C 

It can be seen that the boundary term C just gives a constant contribution to the stress (14). 
The trace of the stress tensor a^k = "Vr + 0" w + cr zz produced by the wedge disclination is 

- = f S{ 1 " r + ^ + /f » M - c }- <15 » 

It is interesting to note that for the case when C = the solution (14) agrees with the result 
obtained by Povstenko [12, 17] in the framework of nonlocal elasticity if we use the identification 
k = 1/{tI). He used the two-dimensional nonlocal kernel which is Green's function of the 
two-dimensional Helmholtz equation. The constant C is determined by the "semi-classical" 
boundary condition on the surface r = R of the cylinder 

a rr {R) = 0, (16) 

which means the absence of external forces on the boundary of the cylinder. So we find for the 
constant 

C = lnR +Y^-^ + Ko(*R) + ^J^( 2 - K 2 R 2 K 2 (KR)y (17) 
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In the limit n^ 1 — > 0, the stress function (12) with (17) agrees with the stress function given 
by [18]. The constant (17) diverges in the limit R — > oo. Therefore, we consider a cylinder 
of finite radius. Note that the cylinder size R is the characteristic screening parameter which 
appears under the logarithm. If we use the limiting expressions for r — > 0, 

1 2 

{nr) 2 

where 7 denotes the Euler constant, we obtain 

MO) ^ <W0) = I o zM = a kk (0) = { In f + 7 - 5 + K (kR) 

' '2- K 2 R 2 K 2 (nR))}. (19) 



#o(«r)->- 7 + lny , X 2 ( K r)^-i + 7 -^, (18) 



Consequently, the stress is finite at the disclination line in contrast to the unphysical stress 
singularity in "classical" disclination theory. The graphs of the components of the stress calcu- 
lated from (14) are plotted over nr in Fig. 1 (with the radius of the cylinder R = 10/ k). With 
R = 10/ k we obtain for (19) the value 

^(0) = , w (0) = ±*„{0) = CTfefe(0 ) * -1.707 ^j^-. (20) 

At the outer boundary i? the non- vanishing components of the stress tensor (14) are 

<W*) = \ °M = Tj^y - ^r^y I 1 - A 2 + 2K ^ R) \ (21) 

With R = 10/ k we obtain for (21) the value (see Fig. 1) 

o w (R) = X -o zz {R) = ^i-j CTfefc (i?) ~ 0.960 (f" - . (22) 

By means of Eqs. (3) and (12) we obtain the modified stress of the wedge disclination in 
the cylinder with radius R in Cartesian coordinates 



>xy 



s^j^-^(>-^.(-))}. 



'"-^{ hp + 5(i^ + Ai ><"'>- }- (23 » 

If we identify k= 1/ yfc (c is the gradient coefficient used by Gutkin and Aifantis) and put C = 0, 
the elastic stress is in agreement with the stress field given by Gutkin and Aifantis [14, 15] in 
the framework of strain gradient elasticity by using the Fourier transform method. 
Using (4) and (12) we find for the elastic strain of the straight wedge disclination 



E xx = i ^ 7 y{(l - 2„)(lnr + *o(«r) - C) + £ + If ^/ 1 (2 - «W a («r)) }, 



K — 



4tt(1 - i/) 
SI 
4tt(1 



{(1 - 2„)(lnr + K (nr) - C) + £ - ^^(2 - «Vx 2 m) } 



£ '" ^ -T^T) ^ 2 ^(-))}" (24) 
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Figure 1: The stress components for a wedge disclination in a cylinder with radius R = 10/ k: 
(a) <r rr , (b) a vv , are given in units of ^f2/[27r(l — v)] and (c) a zz is given in units of ^VIv/[tt(\ — 
v)]. The dashed curves represent the classical stress components. 
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In the case of C = the strain (24) coincides with the result given by Gutkin and Aifantis [14, 
15, 19]. The strain reads in cylindrical coordinates 



E, = - 2*)(lnr + K Q {nr) - C) + -1^2 - K W 2 („r)) J, 



4tt(] 
47r(T" 



E, , = „, J - 2i/)(lnr + Jf («r) - C) + 1 - ^(2 - ^Vi^r)) J. (25) 



The dilatation reads 

^ fe = 27 r(l 1 -^) { (1 ~ 2l/ )( lnr + ^°( Kr ) - C ) + I}' ( 26 ) 

By the help of (18) the expressions of the strain (25) and the dilatation (26) read for the limit 
r — > 0, respectively, 



E rr {Q)=E^{Q) = -^— ) {( \ - >2iA 



, kR 1 » , 

In— +7- 



+ 

and 



^2-(2-k 2 ^2M))}, (27) 



E kk (0) = -g^M { m f + 7 - i + *„(«*) + ^ (2 - }. 



(28) 



The strain and the dilatation at the cylindrical surface r = R are given by 

E VV (R) = £ rr (i?) = -^—^ £ fefe (#) = ^- (l - -i-j + 2K 2 (kR) }. (2!)) 



(1-2!/) ' 4tt 1 K 2 i? 2 



(«i?)j. 



In the conventional disclination theory the torsion tensor (linear version of Cartan's torsion) 
is defined by 

otij = e jki (d k Pii + eum<Pmk) ■ ( 30 ) 

It is the dislocation density in the theory of disclinations (see, e.g., [2,3,20]). On the other 
hand, Anthony [1] called it the disclination torsion. (3ij is the elastic distortion. The ip*j was 
introduced by Mura [20] who called it "plastic rotation" and deWit [2, 3, 21] called this quantity 
"disclination loop density". Using the elastic bend-twist tensor (see, e.g., [22]) 

kij = djUi - if* j: (31) 

with the rotation vector 

= — ^ e ijkfljk, (32) 

Eq. (30) can be rewritten according to (see also [1,21,22]) 

otij = £jki(dkEu + eimikmh) — €jkid k Eu + Sijku — kji. (33) 

Consequently, the elastic bend-twist may be determined from the condition that the dislocation 
density (disclination torsion) has to be zero for a straight wedge disclination 

a xz = ~~2j x ~ dyAj ~ kzx = °' ayz = ~2^~ dxA f ~ kzy = °- ( 34 ) 
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Figure 2: Elastic bend- twist k zlp of a wedge disclination (solid) is given in units of Q/{2ir}. The 
dashed curve represents the classical solution. 



So we find for the elastic bend-twist 

k zx = ^{l - KT\K"i(Kr)j, k zy = ^- ^{l - nrK^nr)}, (35) 



and in cylindrical coordinates 



— - |l - nrKi (nr) j , k zr = 0. (36) 

Obviously, the components of the elastic bend-twist (35) and (36) fulfil an inhomogeneous 
Helmholtz equation 

(l-«- 2 A)% =kij, (37) 

o 

where kij denotes the classical elastic bend-twist. The elastic bend-twist k zip is plotted in Fig. 2. 
It has no singularity at the disclination line. It is k zv (r = 0) = 0. The elastic bend-twist has a 
maximum of k™^ ~ 0.399k£1/[2it] at r ~ 1.1/ n. The boundary condition k zr (R) = is trivially 
satisfied. 

The elastic bend- twist tensor can be decomposed according to (31) into a gradient of the 
rotation vector and an incompatible part. We identify the incompatible part with the disclina- 
tion loop density. It is analogous to the decomposition of the elastic distortion of a dislocation 
into a gradient of the displacement vector and an incompatible distortion (see [6-8]). The 
disclination loop density turns out to be 

<P ' X = £ k2xK ^ kt ) {v>~\ si § n (y)) ' 

fly = ^^K 2 yK (Kr)(^p- | sign(y)) + nS(y) (l - sign(z) [l - nrK^nrj] ) }. (38) 

It contains the angle ip and its form is analogous to the plastic distortion of a screw dis- 
location (see [8]). Here we use a single- valued discontinuous form of ip (see [3]). Only the 
component (p* zy has a <5-singularity at y = like the disclination loop density [3,20] ip* z — 
(0/2) 5(y)(l — sign(x)). The rotation vector reads 

u z = — - nrKiinr)) + ^ sign(y) KrK^nr)}. (39) 
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Figure 3: Rotation vector of a wedge disclination uj z (x,y 
represents the classical solution. 



+0) (solid). The dashed curve 



The far field of the rotation vector (39) agrees with de Wit's expression given in [3]. When 
y — > +0, the expression (39) is plotted in Fig. 3. One can see that the Bessel function terms 
which appear in (39) lead to the symmetric smoothing of the rotation vector profile, in contrast 
to the abrupt jump occurring in the classical solution. It is interesting to note that the size 
of such a transition zone is approximately 12/k which gives the value Q/k for the radius of 
the disclination core. Replacing by 6 in Eq. (39) it can be seen that the rotation of the 
wedge disclination to z has the same form as the displacement of a screw dislocation u z given 
in [8]. Consequently, the elastic rotation (antisymmetric part of the clastic distortion) of a 
wedge disclination is discontinuous due to ip in contrast to the dislocation case where the 
corresponding expressions are continuous (see also [3]). 
Finally, we find for the elastic distortion 



fixx 
fixy 

Pyx 

Pyy 



Ml _ v) { (1 - 2,)(lnr + K ( K r) -C)+^+ ( 2 - «W a («r)) } 

- ^ - nrKi(nr)) + ^ sign(y) KT\Ki(«;r)}, 



+ 



— - nrKi(nr)) + ^ sign(y) KT\Ki(«r)}, 



4tt(1 - v) 



2 / 2 2\ 

(1 - 2i/)(lnr + Koinr) - C) + ^ - - ~\ ' h - k 2 t 2 K 2 {ki 



(40) 



We note that the elastic distortion (40) contains the angle ip in contrast to the dislocation case. 
The elastic distortion satisfies the following relations that the Burgers vector vanishes for a 
wedge disclination 



* {f3 xx dx + f3 X ydy) =0, j> ((3 yx dx + (3 yy dy) = 0. 



(41) 



Here 7 denotes the Burgers circuit. 
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We obtain for the effective Frank vector of the wedge disclination 

fi(r) = <j> (k zx dx + k zy dy) = Q jl - kt\Ki (/cr) j . 



(42) 



This expression is in agreement with the one obtained in [23]. The form of the effective Frank 
vector is similar to the effective Burgers vector of screw and edge dislocations (see, e.g., [6-9]). 
It differs appreciably from the constant value 12 in the region from r = up to r ~ 6/k. In fact, 
we find f2(0) = and Q(oo) = ft. Thus, it is suggestive to take r c ~ 6/k as the core radius of 
the disclination. 

The so-called disclination density tensor of a discrete disclination is defined by [1-3,20-22] 



-tjmndmtPin- 



(43) 



Here the index i indicates the direction of the Frank vector, j the disclination line direction. 
Thus, the diagonal components of 6^ represent wedge disclinations, the off-diagonal com- 
ponents twist disclinations. The disclination density tensor satisfies the continuity condition 
djQij — which implies that disclinations do not end inside the body. From (35) and (38) we 
find for the non- vanishing component of the disclination density of a wedge disclination 



<dzz = K o(Kr). 

Z7T 



(44) 



In the limit as k _1 — > 0, the result (44) converts to the classical expression Q zz — QS(r). It is 
interesting to note that (44) coincides with the expression calculated in [23] . Additionally, the 
disclination density (44) agrees with Eringen's two-dimensional nonlocal kernel used in [10-12]. 

We use the decomposition of the elastic distortion (40) into a compatible and an incompatible 
distortion 



In this way we restore a modified displacement vector (see also 
u x = ~ KrK 1 (nr)) + ^ sign(j/) nrK^nr) 

(1 - 2u) ( lnr - 1 + K {nr) - C) - — 



(45) 



n 

27 



2(1-1/) 



Lp(l — KrKi(nr)) + — sign(y) KrK\{Kr) 



+ 



y 



2(1-1/) 
and the incompatible part 

Pxx = ^ ^KrKi(nr) + K 2 xyK (Kr)(v - 

Pxy = ^K 2 y 2 K (nr)(tp- | sign(y)) , 

Pyx = -^n 2 x 2 K a (nr)(tp - ^sign(y)), 

P ^ 



(1 - 2v)(lnr - 1 + Ko(nr) - C) 



(2- n 2 r 2 K 2 {nr)) 



(2- K 2 r 2 K 2 {nr)) 



(46) 



|sign(y))|, 



(47) 



" ~ 2~7r { KrKl ( Kr ) _ K 2 xyK Q (nr)(tp - ^ sign(y)) - TrS(y)x(l - sign(x) [l - ktK x (/cr)] ) | . 



Eq. (47) fulfils atij = ejkddkPii + etimPmk) = and Eqs. (46) and (47) satisfy u = d[ y u x ] +P[ xy \- 
The 5-term in (47) has a similar form like the plastic strain or dislocation loop density of a 
wedge disclination [3, 20, 21] E y = (5* = -(0/2) 8(y)x(l - sign(x)). 
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Using the gauge theory of defect [24], which is an 7S'0(3)-gauge theory (I SO (3) = T(3)2> 
5*0(3)), we are able to decompose the incompatible distortion (47). Namely, the incompatible 
distortion takes the form [24, 25] 

Pij = <Pij + £ikiW kj xi, (48) 

where faj and Wij are the translational and rotational gauge fields, respectively. We obtain 
for (47) the following decomposition 

Pxx = 4>xx - yW zx , 

Pxy = ~yW zy , 

f} yx = xW ZXl 

(3yy = 4>yy + xW z y , (49) 

into the translational gauge field 

4>xx = <i>yy = i^ KrKi(nr), (50) 

and the rotational gauge field 

W ZX = ~<P* ZX , W Z y EE ~ V * zy . (51) 

Therefore, the negative disclination loop density (38) is equivalent to the rotational gauge 
potential (51). 

In conclusion, the field theory of elastoplasticity has been employed on the consideration 
of a straight wedge disclination. We were able to calculate the elastic and plastic fields. We 
found that the elastic stress, elastic strain, elastic bend-twist and disclination density are con- 
tinuous and the displacement, plastic distortion, rotation and the disclination loop density of 
the wedge disclination are discontinuous fields. Exact analytical solutions for all characteristic 
field quantities of a wedge disclination in a cylinder have been reported which demonstrate 
the elimination of "classical" logarithmic singularities at the disclination line. All logarithmic 
terms are influenced by the "semi-classical" boundary term C (17). In addition, the disclina- 
tion core appears naturally as a result of the smoothing of the rotation vector profile. For an 
infinitely extended body (C = 0) the elastic stress of a wedge disclination calculated in the 
field theory of elastoplasticity agrees with the stress calculated within the theory of nonlocal 
elasticity and strain gradient elasticity. The reason is that in all three theories the fundamen- 
tal equation for the elastic stress has the form of an inhomogencous Hclmholtz equation (see 
Eq. (1)). The boundary- value problem of a wedge disclination in a cylinder considered in this 
paper should be help in studies of mechanical behaviour of nano-objects including nanotubes 
and nanomembranes and of disclinated nanoparticles of cylindrical shape (nanowires). Finally, 
we note that one observes an interesting relation between the wedge disclination and the screw 
dislocation. Namely the rotation (39), elastic bend-twist (36), effective Frank vector (42) and 
the disclination density (44) of a wedge disclination have the same form as the displacement, 
elastic distortion, effective Burgers vector and the dislocation density of a screw dislocation 
given in [7, 8] when the Frank vector is replaced by the Burgers vector. 
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